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We investigate the entanglement between a spin and its environment in impurity systems which 
exhibit a second-order quantum phase transition. As an application, we employ the spin-boson 
model, describing a two-level system (spin) coupled to a subohmic bosonic bath with power-law 
spectral density, oc ui" and < s < 1. Combining Wilson's Numerical Renormalization Group 

method and hyperscaling relations, we demonstrate that the entanglement between the spin and 
its environment is always enhanced at the quantum phase transition resulting in a visible cusp 
(maximum) in the entropy of entanglement. We formulate a correspondence between criticality and 
impurity entanglement entropy, and the relevance of these ideas to Nano-systems is outlined. 

PACS numbers: 73.43.Nq, 72.15.Qm, 03.65.Ud 



Quantum mechanical systems can undergo zero- 
temperature phase transitions upon variation of a non- 
thermal control parameter, where the order is destroyed 
solely by quantum fluctuationsii In this Letter, we re- 
explore second-order quantum phase transitions in impu- 
rity models which involve a spin coupled to a dissipative 
environment and which display both a localized and a de- 
localized phase for the spin.— To better characterize those 
quantum phase transitions, we are prompted to exam- 
ine the entropy of entanglement shared between the spin 
and its environment,-'- As a concrete example, we em- 
ploy the subohmic spin-boson modelS. describing a two- 
level system with environmental dissipation (stemming 
from a lossy RLC transmission lineS. or from 1// noise*') 
and which allows a direct measurement of the entropy 
of entanglement. We demonstrate that the ground state 
is strongly entangled at the quantum phase transition. 
The analysis of the enhancement of entanglement near a 
quantum critical point is of great current interest.— i^i^ 

More precisely, we show that those second-order impu- 
rity quantum phase transitions are always accompanied 
by a maximum (cusp) in the entropy of entanglement 
and that the latter exhibits universal scalings. The en- 
tanglement entropy will also allow us to establish impor- 
tant connections between impurity entanglement, quan- 
tum decoherence (or strong reduction of the quantum su- 
perposition of the two spin states) when approaching the 
phase transition from the delocalized phase, and rapid 
disentanglement in the localized or classical phase for the 
spin (the spin is rapidly frozen in one classical state due 
to dissipation). Our phase diagram is shown in Fig. 1. 

The simplest example of a quantum impurity coupling 
to a bath is the paradigmatic spin-boson models^ 



J^SB — ~-^(^x + t;'^z 
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^A„(a„ + at), (1) 



where cr^ and cr^ are Pauli matrices and A is the tunnel- 
ing amplitude between the states with cr^ = ±1. In the 
following, we will assume that A > such that (fix) > 0. 
Here, Hose is the Hamiltonian of an infinite number of 
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FIG. 1: (color online) Unification of entanglement and criti- 
cality at a second-order impurity quantum phase transition. 



harmonic oscillators with frequencies {wn}, which couple 
to the spin degree of freedom via the coupling constants 
{A„}. The heat bath is characterized by its spectral func- 
tion, J{iLj) = TrJ^n ^nH'^n — uj) = 2tt aoj].^ uj'' , wherc Wc 
represents a cutoff energy and the dimensionless param- 
eter a measures the strength of the dissipation. 

The essential physics contained in this model is the 
competition between the amplitude for tunneling be- 
tween the two spin states (leading to a "delocalized" 
phase) and the effect of the bath which tends to "lo- 
calize" the system in one or other of the spin states. 

The special value s = 1 represents the case of ohmic 
dissipation, where the analogy with the Kondo model 
applies and a quantum Kosterlitz-Thouless (KT) transi- 
tion separates the delocalized phase at small a and the 
localized phase at large ai^ The phase transition occurs 
at Q!c = 1-1- ©(A/wc)'^ The delocalized region corre- 
sponds to the antiferromagnetic Kondo model with a 
Fermi-liquid ground state while the localized region cor- 
responds to the ferromagnetic Kondo model where the 
spin is frozen in time. The spin magnetization (cr^) jumps 
by a non- universal amount —1-1- ©(A/wc) at the KT 
transition (for e = 0+).^^ The entropy of entanglement 
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also contains valuable information; it shows a plateau 
at maximal entanglement for a > 1/2 in the delocal- 
ized phase and it drops (roughly) to zero at the KT 
transitioufi^ii^ii^ Note that this plateau exemplifies the 
intimate connection between maximal entanglement and 
quantum decoherence since a = 1/2 marks the dynam- 
ical crossover from damped oscillatory to overdamped 
behavior.^ 

Below, we rather focus on the subohmic spin-boson 
model (0 < s < 1). A quantum critical point Ac (a) 
still separates the localized (A < A^) and the delocal- 
ized phase (A > Ac) where Ac ^ when a ^ 0. A 
second-order phase transition has been well-established 
for all < s < 1 using Wilson's Numerical Renormaliza- 
tion Group (NRG) method^-ii^ and through an analogy 
to classical spin chains^^ (nevertheless, when interactions 
become too long-ranged in time, i.e., for s < 1/2, the 
quantum-classical mapping fails^) . The longitudinal spin 
magnetization (az) is generally used to characterize the 
phases and phase transitions.^" On the other hand, the 
quantum critical point is an "interacting" fixed point for 
all < s < 1 which motivates us to examine the entangle- 
ment between the spin and its environment thoroughly. 
When a bipartite quantum system AB is in a pure state 
lip), the entanglement between subsystems A and B is 
unambiguously described by the von Neumann entropy 
or entanglement entropy E which is calculated from the 
reduced density matrix pA or pB, Pa(b) = TrB(A)|'0)(V'|'^ 



E = -Tr {pA log2 pa) = -Tr {ps logs pb) ■ 



(2) 



When either subsystem A or B is a spin-i sys- 
tem, the entropy of entanglement E can be rewritten 
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E 



-P+ log2 P+ - P- log2 p- , where p± 



as," 

^1 ± {cTx)^ + (o'z)^j /2; note that {ay) = because the 
Hamiltonian Hsb is invariant under ay — > —ay. 

Since exact analytical methods (such as Bethe Ansatz) 
are restricted to the ohmic caser^^ we employ the 
bosonic NRG to compute the entanglement entropy be- 
tween the spin and its environment. We follow the same 
procedure as in our Ref. 16. In order to ensure the con- 
vergence of the results in the localized phase, we have 
used the NRG parameters A = 2 (logarithmic discretiza- 
tion), Ns — 150 (lowest energy levels kept), and Nf, = 8 
(boson states; except the 0th site for which A^f,o = 500). 
The results converge for N !v 30 sites (and we keep until 
40 sites). In NRG calculations, we Rx ojc — ^ and A is 
normalized to loc. We will also apply the following scaling 
ansatz for the impurity part of the free energy. 



F = ^/(|A-Ac|/Tl/^e^-^) 



(3) 



to relate the critical exponents associated with E to other 
critical exponents such as the correlation length exponent 
u. Even though the temperature T is introduced for the 
scaling analysis the entanglement entropy E (which is 
defined for a pure state) is evaluated at zero temperature. 
The crossover from the quantum critical regime to one or 
other of the stable regimes, defines the energy scale e* in 
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FIG. 2: (color online) NRG results for the entanglement at 
T = e = in the subohmic case s = 0.5, which can be realized 
through a qubit coupled to an RLC transmission line; the 
arrow marks the position of the quantum phase transition. 
The relevant susceptibilies Xz and X-L are presented in insets. 

Fig. 1 that vanishes at Ac as e* oc |Ac — Aj*"^, and for the 
spin-boson model we obtain b ~ {I + s)/2. The ansatz 
([3]) is well justified when the fixed point is interacting; — 
for a Gaussian fixed point the scaling function would also 
depend upon dangerously irrelevant variables. 

Note that the entanglement entropy E is different from 
the impurity entropy Simp which is rather evaluated as 
entropy of the system with impurity minus entropy of 
the bath alone. The delocahzed phase obeys Simp{T — > 
0) = since the ground state is non-degenerate while the 
localized phase yields Simp{T — 0) = ln2. On the other 
hand, E goes rapidly to zero in the localized phase (prod- 
uct state) and is finite in the delocalized phase (entangled 
state); however, i? — > when A 3> Wc- Our NRG data of 
Figs. 2 and 3 show that E exhibits a cusp-like behavior 
(maximum) at the phase transition. This is distinguish- 
able from the ohmic case where \{az)\ = 1 — ©(A/tJc) at 
the KT transition which irrefutably results in E Oj ^'^i^^ 

In fact, there is a simple way to conceive that the 
maximum of E coincides with the quantum phase transi- 
tion for the subohmic situation. In the delocalized phase, 
since the longitudinal magnetization {az) = at e = 0: 
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This expression is valid at Ac. Since X-L = d{o'x)/dA 
and {ax) are positive quantities, dE/dA < in the delo- 
calized phase. In the localized phase, E is controlled by 
the finite longitudinal magnetization and by the suscep- 
tibility Xz = -d\{az)\/dA > (see inset in Fig. 2)-^ 
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Thus, dE/dA > in the localized phase. Eqs. (6) and 
(7) imply that the entanglement entropy is maximum at 
the phase transition. Zero-temperature impurity criti- 
cal points can show a residual "fractional" entanglement 
entropy (which depends on the dissipation strength a). 
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FIG. 3: (color online) The entropy of entanglement E, ob- 
tained from NRG, displays a cusp-like behavior for all < 
s < 1. For s — > 1, _E becomes rapidly suppressed at A~ 
which is a reminiscence of the ohmic case (KT transition). 

Let us emphasize at this point that the case s — 1/2 
is of particular interest since it can be realized through a 
charge qubit (dot) subject to the electromagnetic noise of 
an RLC transmission line^ the mapping onto the spin- 
boson model with s = 1/2 assumes frequency to <^ R/L, 
whereas an LC transmission line would mimic the ohmic 
casei^i^'^ The entanglement entropy can be accessed 
through charge and persistent current measurements, 
corresponding to (cr^) and (ax) respectively.^^ Moreover, 
A represents the tunneling amplitude between the dot 
and the lead or the Josephson energy of the junction and 
thus is a tunable parameter. A gate can be used to con- 
trol e. The noisy charge qubit is an ideal Nano-candidate 
to demonstrate the existence of entanglement 

Now, we seek to demonstrate that E exhibits universal 
scaling even though the entanglement is two-sided, so 
that two numbers are necessary to specify E (one for 
each of the two ways of approaching Ac). Near Ac, the 
transverse spin susceptibility x± — 29^F/(9A)^ obeys 

X±(A)-X±(Ac) + c+/_|A-Ac|«, (6) 

and from the free energy defined in Eq. ([3]), C = ^ 2. 
For the subohmic spin-boson model, one finds ly > 2 for 
all < s < 1, ensuring that x± does not diverge at the 
transition. Taking into account that (ax) is continuous 
at the transition, Eq. ([6]) thus implies that E always rises 
linearly for A A+ (A^ means that we approach the 
quantum critical point from the delocalized (localized) 
region), as confirmed by Fig. 3. Note that the coeffi- 
cients c+ and c_ can be different in the delocalized and 
in the localized phase and this explains for example the 
jump occurring for s = 1/2 where v = 2 (inset in Fig. 
2). Through the NRG, we also check that c+ < 0, em- 
phasizing that in the delocalized phase x± substantially 
increases at A+, and that x± shows a clear singularity 
for 2 < 1/ < 3 or 1/3 < s < 0.94. For all < s < 1, this 
strongly underlines the duality between the enhancement 



FIG. 4: (color online) -E(A) at s = 0.5 and a = 0.1 for several 
values of e. Inset: universal scaling in the critical region for 
three different values of Ac (from NRG we obtain k ~ 0.6 and 
the exact exponent is k = 2/3 = 0.66...). 

of entanglement and the strong reduction of the two-spin 
state quantum superposition near the phase transition. 
In the localized phase, we obtain the scaling behavior: 

X.(A)cx|A-Acri+'^(i-^)/2 + a; (7) 

here a ^ when 1/(1 — s)/2 > 1, and we identify a = 
Xz(Aj). For A > Ac, Xz — 0. The correlation length 
exponent can be obtained analytically for s and 
s — > 1 through Renormalizations Group expansions. The 
analogy with classical spin chains for s — > 1 leads to^i 
— \/2{l — s) whereas at small s, one finds^ l/v ~ s. 
Here, Xz diverges at A~ for s > 1/3 and /3 = i/(l — s)/2 is 
the critical exponent associated to {(Jz)- For A < Ac and 
s > 1/3, from Eqs. (7) and (9), we find the relation^^ 

£(Ac)-S(A)oc|A-Acr(i-^^ (8) 

The decay of the von Neumann entropy E in the localized 
phase is faster than linear for all s > 1/2 (Fig. 3) and 
the behavior becomes strictly linear at s = 1/2, as shown 
in Fig. 2. It is certainly relevant to notice the parallel 
between impurity entanglement in a dissipative environ- 
ment and single-site entanglement in quantum critical 
spin chains such as the anisotropic XY chain.-i Now, we 
shall discuss the scaling of E with the longitudinal field. 

Integrating out the boson degrees of freedom induces 
a long-range interaction in time which results in the fol- 
lowing term in the action;^ Si^t = / drdr' a2{T)g{T — 
t')(Tz{t'), where f/(r) oc 1/t^+'' at long times. Assuming 
that the dynamics of cr^ at the critical point is essentially 
determined by Smt and by the local field e we then derive 
{az){e,Ac) oc |e|i/'', with the exponent (5 = {l+s)/{l-s). 
This is consistent with our NRG results which predict 
that the local susceptibility Xz = d\{az)\/de at the quan- 
tum critical point diverges as T^'*. In the small s limit, 
this leads to 1/S ~ 1 — 2s -I- 0{s'^); this can be recov- 
ered by resorting to a small s expansion.^ Using (ax) = 
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FIG. 5: (color online) Exponent n versus s from NRG. For all 
< s < 1, K = 2/(5 and for < s < 1/2, k = 2/5 = 1/5. 



-2dF/d_A, we obtain K)(e,A,) - K)(A,) oc -|e|i/^ 
with 1/(5 = — l/i^)- To maintain consistency with 

our notations, we identify {(Jx){£ — 0, Ac) = {(Tx){Ac), 
and similarly for other quantities (observables). 

When s is close to one, the critical point is defined by 
a — I and A^ = 1 — s. This ensures that (ax) is small 
at the transition and evolves slowly with e; at s = 1, 
around the phase transition, one gets the exact expres- 
sion {(Tx){e = 0) = A/[u!c{2a—l)]}'^ Thus, when s is close 
to one, the dependence of i? on e mainly stems from ((Tz): 

i?(e,Ac)-i?(A,)cx-|er, (9) 

and K = 2/S = 2(1 - s)/(l + s). In fact, since l/S{s) > 



2/6{s), this scaling relation remains valid for all < s < 
1, as shown in Figs. 4 (inset) and 5; one can always 
expand p± (e) — p±{e = 0) i: me^/* at small e and m > 
to satisfy deE{e) < (the field e favors a product state). 
On the other hand, since at small s the critical exponent 
1/ obeys = 1/s, one also gets 1/6 — 2/6 — k, which is 
well verified through the NRG for < s < 1/2 (Fig. 5). 

For A > Ac, the NRG results predict {a^){e,A) oc e 
and {ax){e, A) — {ax){A) oc —e^. Thus, E decreases as 
similar to the ohmic caseJ^ Since 2/6 < 2 for all < s < 
1, this implies that for a given e 7^ 0, the maximum of 
entanglement occurs at the value of A > Ac which lies in 
the crossover between the delocalizcd and the quantum 
critical regime (Figs. 1 and 4). The dclocalized phase is 
quite robust to the application of a field e. For A < Ac, 
in contrast we find a linear decrease of E with e. 

In conclusion, we have shown that the entanglement 
between a spin and its (bosonic) environment is always 
enhanced at a second-order quantum phase transition. 
The concept of entanglement entropy allows us to estab- 
lish important connections between impurity entangle- 
ment, strong reduction of the quantum superposition of 
the two spin states when approaching the phase transi- 
tion from the delocalized phase, rapid disentanglement 
in the localized phase, and criticality. Our theoretical re- 
sults can be tested experimentally through a charge qubit 
coupled to a lossy RLC transmission line. These results 
may also be relevant for heavy fermion systems which 
might develop a similar "local" criticality.^^ 
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